Influence of the trap shape on the detection of the superfluid-Mott transition 
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The coexistence of superfluid and Mott insulator, due to the quadratic confinement potential in 
current optical lattice experiments, makes the accurate detection of the superfluid-Mott transition 
difficult. Studying alternative trapping potentials which are experimentally realizable and have a 
flatter center, we find that the transition can be better resolved, but at the cost of a more difficult 
tuning of the particle filling. When mapping out the phase diagram using local probes and the local 
density approximation we find that the smoother gradient of the parabolic trap is advantageous. 
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Ultracold atomic gases in optical lattices provide a tun- 
able and well-controlled realization of strongly correlated 
quantum lattice models such as the bosonic or fermionic 
Hubbard models [lj]. Unlike strongly correlated materi- 
als, where the microscopic model is often unknown and 
the situation is complicated by additional interactions 
such as phonons and long-range Coulomb interactions, 
cold atoms provide a clean realization of prototypical 
models of strongly correlated systems. 

In a seminal paper, Greiner et al. 01 demonstrated 
the loss of coherence in a Bose-Hubbard model as one 
tunes the interaction parameters from a superfluid Bosc- 
Einstein condensate (BEC) to a bosonic Mott insulator 
(MI) in a parabolic trap. While this experiment beauti- 
fully demonstrates the BEC and MI phases, it remains at 
the qualitative level of demonstrating the existence of the 
two phases. To use ultracold atomic gases as a quantum 
simulator for strongly interacting quantum systems one 
has to go an important step further and quantitatively 
determine the phase diagrams. Apart from devising new 
detection techniques, there are two major challenges to 
be solved: thermometry and the influence of the trap 
shape. Here we focus on the latter aspect. 

While the Greiner et al. experiment identifies the BEC 
and MI phases, the "transition" is a smooth and broad 
crossover between these two limits instead of a true phase 
transition due to the inhomogeneous parabolic trapping 
potential The coherence fraction (i.e., the bright- 
est spot of the interference pattern) does not display 
any distinct features upon emergence of a central Mott 
plateau, but decreases smoothly over a broader range 
than expected for the uniform system. Similarly, the 
FWHM barely changes near the Mott onset, but increases 
abruptly for much deeper lattices. While for the homo- 
geneous lattice the transition can be more or less accu- 
rately determined using the coherence fraction and the 
FWHM, no such approach is possible in the presence of 
a parabolic trap, and the Mott onset cannot simply be 



inferred from time-of-flight experiments 3. 

We can attempt to obtain better information about the 
phase transition of the homogeneous system by study- 
ing the crossover in a trap if we use "flatter" traps with 
steeper walls, mimicking a homogenous system with hard 
walls, since in such a trap the crossover is expected to be 
sharper and closer to the true phase transition. 

In this Letter we investigate the effect of the shape 
of the trapping potential V(r) on the detection of the 
BEC-MI transition in two-dimensional bosonic lattice 
gases, described by the single-band Bose Hubbard Hamil- 
tonian 
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For the hopping amplitude t and the on-site repulsion 
strength U we use the microscopic parameters in the 
standard tight-binding approximation Q for 87 Rb atoms 
with an s-wave scattering length of a s = 102(6)ao [f|, 
where ao is the Bohr radius. We obtain a two- 
dimensional lattice by integrating over the strong har- 
monic confinement in the z-direction with oscillator 
length a z = 53.5nm and assume a lattice laser with wave- 
length A = 820nm, which gives a lattice spacing a = A/2. 
The recoil energy En = h 2 / (2mA 2 ) is the natural unit to 
express laser intensities. 

We will investigate three different trapping potentials: 
parabolic V(r) = V^/a) 2 , quartic V(r) = V±(r/a) A and 
flat potentials with "light sheet" walls. The quartic po- 
tential can be realized either using supergaussian laser 
beams or by compensating the typically quadratic part 
of the trapping potential with a blue-detuned laser. An 
even flatter trap center can be obtained in a "light sheet" 
trap using repulsive walls of blue-detuned Gaussian laser 
beams. The waist of the light sheet laser beam is set to be 
12 lattice spacings, i.e. a = 6a. The trapping potential 
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FIG. 1: (Color online). The three trapping potentials used 
in our comparison along one direction. The scale is different 
for each of the potentials, in particular the walls of the light 
sheet potential are much steeper than suggested in the figure. 



of these light sheets is 
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where r = (x,y) is the distance from the center of the 
trap and we have placed light sheets parallel to the x 
and y-axis centered along x = ±r c and y = ±r c . These 
trapping potentials are schematically shown in Fig. [1] 

In all cases we assume that the whole confinement is 
coming from the external potential only and that there 
is no contribution from the lattice laser waists or the z- 
direction confining potential. This can in principle be 
done by carefully compensating the Gaussian trapping 
potential of the waist of a red-detuned lattice laser with a 
blue-detuned dipole trap. The required homogeneity for 
the lattice potential depends on the energy and length 
scales of the problem. Globally, the lattice potential 
should be flat to much better than the local chemical 
potential at the transition, which for our 87 Rb exam- 
ple is set by U and should be smaller than E-r over the 
50 x 50 sites in the sample. Locally, the site-to-site poten- 
tial should be flat to much better than t at the transition, 
thus variations should be much smaller than 0.05£/a over 
the a = 410 nm lattice spacing. We note that for "super- 
exchange" physics which scales as t 2 /U, these require- 
ments are much more stringent. In any case, using larger 
gaussian beams provides significantly more homogeneous 
lattices, at the expense of increased required laser power. 

We investigate the superfluid-Mott transition in these 
trap ping potentials using Gutzwiller mean-field theory 
[l2j at zero temperature with a fixed particle num- 
ber N. We miminize the ground state energy Eg = 
(G\H\G) of the mean-field Gutzwiller ansatz \G) = 

Hi (j2^=o frf* \ n )ij I where \n) i denotes the Fock state 
with n particle numbers at site i. While Gutzwiller mean- 
field theory does not give the exact value of the critical 




point, a comparison with an accurate phase diagram 13 1 



shows that the mean-field approach captures the physics 
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FIG. 2: (Color online). Condensate fraction (main figure) 
and its second derivative (inset) as a function of U for a num- 
ber of different fillings and a light sheet trap with amplitude 
Vls = 18-E/r. The light sheet is centered at r c — (0, ±38a) 
and (±38o, 0) and there are 2550 bosons. For comparison 
we include the condensate fraction in a homogeneous system. 
The inset shows the second derivative n" = (1/N)d 2 n(k = 
0)/d(U/t) 2 for the same curves. 



of the BEC-MI phase transition and is reliable for the 
present purpose of qualitatively comparing the effects of 
different trapping potentials. 

Estimating the transition from the coherence frac- 
tion - In a homogeneous system the Mott transition 
is clearly observed by the vanishing condensate fraction 
f c = n(k = 0)/N. In a trapped system the smooth 
crossover due to phase coexistence makes the definition of 
the transition point somewhat ambiguous. To estimate 
it, we note that there is a cusp in n(k = 0) at the tran- 
sition point in the homogeneous system. We determine 
the " transition" point in the trapped system as the point 
where the second derivative of the condensate fraction 
reaches a maximum Q. As we can see form the inset in 
Fig. [2] this gives a 'transition point' in good agreement 
with that in a homogeneous lattice. 

We find a strong dependence of the critical point on the 
particle number, as shown in Fig. [3J which is important 
in view of the difficulty of controlling particle number in 
experiments. Using the minimum of the transition U/t as 
a function of particle numbers as an estimate for the tip 
of the Mott lobe, we find only a small dependence on the 
trap shape. The condensate fraction and the estimate for 
the transition point at these minima are shown in Fig. |4] 
Changing a parabolic trap to a quartic or light sheet trap 
decreases the deviation compared to the homogeneous 
system from 1.2% to 0.6% and 0.3%, respectively. 

Looking at density profiles in Fig. [5] we see that when 
a central Mott insulating region is formed at the tip of 
the Mott lobe, the width of the surrounding supcrfluid 
region in the light sheet setup is considerably smaller 
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FIG. 3: (Color online). The estimate for Mott transition as 
a function of the total number of particles N for parabolic, 
quartic and light sheet traps. The values for the trapping po- 
tentials are chosen as V 2 = 2.9 x 10~ 4 S R , V 4 = 3.3 x 10" 7 £ R , 
and Vls = l&En in order to have Mott domains of comparable 
size and the minimum to occur at similar particle numbers. 
The dashed line indicates the critical point as obtained in the 
same mean-field approximation. 



than for the quadratic and quartic trap. The density 
in the supcrfluid region also drops much faster as can be 
expected from the shape of the potential. This is reflected 
in the condensate fraction curves of Fig. [2 where we also 
compare to the homogeneous case. When U/t is large, 
the light sheet potential has a much lower condensate 
fraction than the quartic or parabolic potential, because 
the supcrfluid edges are almost absent. The gradient 
change in n(k = 0)/N with U is more abrupt for the light 
sheet potential than for more curved potentials. Other 
quantities than the coherent fraction that are also related 




FIG. 4: (Color online). Comparison between the optimal 
curves for the parabolic, quartic and light sheet potential. 
It is very hard to determine the transition point with the 
parabolic trap. The light sheet trap gives the best estimate 
for the phase transition. 
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FIG. 5: (Color online). Density profile (left, full line) and 
trapping potential (right, dashed line) for U/t at the tip of 
the Mott lobe for TV = 2550 particles. 



to the TOF interference pattern measurements, such as 
the FWHM of the first peak or the visibility (n max — 
imin)/(fimax + "■min) are also insensitive probes for the 
superfluid-Mott crossover in a parabolic trap. We have 
checked with quantum Monte Carlo that for the light 
sheet potential, we obtain FWHM and visibility curves 
that are closer the homogeneous lattice results - again 
confirming that a flatter trap is advantageous if atom 
number and temperature can be well controlled. 

The main advantage of the light sheet trap is thus the 
much sharper signal (compare the bends in the curves of 
Fig. HJ. The light sheet trap is, however, also the most 
sensitive one to particle number fluctuations. While the 
transition in the parabolic trap varies by less than 1% 
when changing N by 10%, the transition in the light sheet 
trap varies as much when changing N by only 1%. It is 
clear that quasi-flat traps such as the light sheet trap are 
extremely sensitive to the particle number; the transition 
to the Mott phase is largely density-driven; a big value of 
U/t alone is not sufficient for observing the Mott state in 
flat traps. Another challenge comes from entropy consid- 
erations, which must be taken into account since ramping 
up the optical lattice is to a good approximation an adi- 
abatic and not an isothermal process. Since for the com- 
mensurate, homogeneous lattice the Mott gap opens as 
~ U ft, temperature has to increase proportionally. For 
the parabolically trapped Bose-Hubbard model, temper- 
ature almost scales with U/t, but with a pre- factor that 
is typically an order of magnitude lower thanks to the 
non-zero superfluid volume fraction flij . In the latter 
case, temperature effects play only a minimal role, but 
for the light sheet potential they might be important and 
require a colder initial state that has a much lower en- 
tropy. However, once 1% accuracy in atom number can 
be achieved experimentally, a quartic or light sheet trap 
will give a substantially easier to observe experimental 
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FIG. 6: (Color online). Phase diagram obtained using the lo- 
cal density approximation for the parabolically trapped two- 
dimensional Bose-Hubbard model. The phase boundary is 
determined by looking where the density deviates from unity 
by a threshold value e. We plot the results obtained for var- 
ious values of e: e = 0.001 to e = 0.010. The procedure fails 
close to the tip of the lobe. 



signature and a better estimate for the critical point. 

Obtaining the phase diagram from the local density - 
A new and promising detection technique is the high- 
resolution quantum gas microscope [l5| , giving access to 
in situ measurements. Using the local chemical potential 
approximation, /j, Q s(i) = y,—V(i), where \i is the chemical 
potential, one can map out the grand canonical phase 
diagram since for every value of U/t one obtains a slice 
of the (U/t, n/U) phase diagram. We want to investigate 
which trapping potential is best suited for obtaining the 
phase diagram this way. 

Again the inhomogeneous trap smears out the transi- 
tion and the local chemical potential approximation is 
not valid near the phase boundaries There is thus a 
degree of arbitrariness (or uncertainty) in the procedure: 
looking at the density profiles of Fig. [5) we have to define 
where the Mott insulator ends. 

One criterion could be the deviation of the density from 
an integral value. Choosing a threshold — 1 < e as 
the criterion for detecting the phase transition to the 
n = 1 Mott lobe, we plot the obtained phase diagram 
in Fig. [5] Plotting the estimates for various values of 
£ in a parabolic trap in Fig. [5] gives us an idea of the 
uncertainty of the phase boundary. The resulting phase 
diagram is in reasonably good agreement with the homo- 
geneous model. Near the tip of the lobe our procedure 
does not work: we do not obtain any data points there 
because the density goes initially down along a line of 
constant chemical potential starting at the critical point. 
It is worth noting that had we determined the phase dia- 
gram based on density variations the tip of the lobe could 
be described accurately. 

It turns out that the parabolic trap is best suited for 
obtaining the phase diagram using local probes since the 
trapping potential increases the slowest with distance. 
This gives the smallest errors for the local potential ap- 



proximation and provides a finely spaced grid of chemical 
potentials over a wide region. A quartic or light sheet 
trap results in substantially larger deviations from the 
homogeneous result, since using local probes the phase 
diagram is determined from the edge of the Mott region, 
which here is in the steep boundary zone of the trap and 
not in the flat bottom. 

Conclusions - We have studied the two-dimensional 
Bose-Hubbard model with a parabolic, quartic and light 
sheet trapping potential. We have seen that the light 
sheet potential trap gives a much clearer signal for the 
superfluid-Mott crossover. It is however much more sen- 
sitive to the total number of atoms in the system since 
commcnsurability is important when the inner part of 
the trap is almost fiat. Nevertheless, our values for (U/t) 
at the crossover are in good agreement with mean-field 
theory for a homogeneous lattice. 

Measuring the local density combined with a local 
chemical potential approximation allows an easier deter- 
mination of the phase diagram. This procedure allows 
mapping out the full phase diagram in the (U/t, n/U)- 
plane and a parabolic trap is preferred over flatter traps 
since the flatter slope gives access to more values of the 
chemical potential and improves the accuracy of the local 
chemical potential approximation. The obtained phase 
diagram is in good agreement with the homogeneous 
phase diagram. 

We thus conclude that currently the most promising 
procedure to quantitatively determine the phase diagram 
of the Bose-Hubbard model are local measurements in 
a harmonic trap - it will still be a challenge to obtain 
high accuracy on the phase boundary because of errors 
due to the local density approximation. Quantum Monte 
Carlo (QMC) simulations of three-dimensional systems 
with realistic parameters are in progress to estimate the 
accuracy of this approach. The position of the tip of the 
Mott lobe could be determined more accurately using a 
global measurement of n(k), but the required accuracy 
in particle filling and in measuring the momentum dis- 
tribution to find the largest curvature is a challenge for 
the experiments. 
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